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Gravity-mode pulsations in subdwarf B stars: a critical 
test of stellar opacity 
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ABSTRACT 

The identification of non-radial g-mode oscillations as the cause of variability in cool 
subdwarf B stars (PG1716 variables) has been frustrated by a 5 000 K discrepancy be- 
tween the observed and theoretical blue edge of the instability domain ({Fontaine et alJ 
2003). A major component in the solution to this problem has been identified by (a) 
using updated OP instead of OPAL opacities and (b) considering an enhancement of 
nickel, in addition to that of iron, in the driving zone. The reason for this success is 
that, in OP, the "Fe-bump" contributions from iron and nickel occur at higher tem- 
peratures than in OPAL. As well as pointing to a solution of an important problem 
in stellar pulsation theory, this result provides a critical test for stellar opacities and 
the atomic physics used to compute them. 

Key words: atomic processes, radiative transfer, stars: interiors, stars: oscillations, 
stars: horizontal branch, stars: early-type 



1 INTRODUCTION 



Ijefferv fc Said ll2006l) recently reviewed the excitation of 
pulsations in certain low-mass stars by so-called "Fe- 
bump instability". This instability is caused by a signifi- 
, cant contribution to stellar opacity from M-shell electrons 
in iron-group elements at temperatures around 200 000 K 
jRoeers fc Iglesias! Il99l iThe Opacity Project Teamlll995D . 
It provides the driving mec hanism for radial pulsations in 
certain extreme helium stars iSaioll993f) and p-mode oscilla- 
tions in hot subdwarf B stars (EC 140 26 variables: after pro- 
totype EC 14026-2647 = V 361 Hya) llCharpinet et al.lll996t 

ll997tlKilkennv et alJll997h . 

The Fe-bump opacity mechanism is effective in these 
stars because of an increase in the contrast between opacity 
due to iron-group elements and opacity from other sources. 
In the case of extreme helium stars, the background hy- 
drogen opacity is sup pressed. In the case of EC 14026 stars, 

Ijefferv fc Said ( 2006) demonstrated that, for radial and non- 
radial p-mode oscillations, the observed boundaries of the 
instability strip can be explained only by increasing the 
iron abundance (Ap e ) without increasing the heavy element 
abundance (Z) as a whole. It was already well known that 
radiative forces act differentially on the ions in the stellar 
envelope such that substantial chemical gradients are es- 
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tablis hed over a diffusion time scale ~ 10 5 y llMichaud et alJ 
1989). The consequent levita tion and accumulation of iron in 
layer s at around 200 000 K dChaver. Fontaine fc Wesemaell 
Il995l) enhances the Fe-bump opacity sufficiently to excite 
pulsations in ab out 10% of sdB stars within the EC 14026 
instability zone llCharpinet et al.ll200lT) . 

The discovery of oscillations with periods of a few 
hours in many cool sdB stars (PG1716 variables: after 
PG 1716+4 26) has presented a challenge to stellar pulsa- 
tion theory jGreen et alJ2005t . While the radiatively-driven 
diffusion of iron can still operate in these stars, p-modes 
wer e reported to be stable in the chemically stratified mod- 
els llCharpinet et all 1200 ll) . On the other hand, non-radial 
g-modes of high radial order (fc > 10) and high spherical 
degree (I > 3) were found to be unstable, but only in mod- 
els of sta rs cooler than those in which variability has been 
detected llFontaine et alj|2003h . 

Ijefferv fc Said (120061) found that, for homogeneous 
models of blue horizontal branch stars, non-radial g-modes 
of high radial order (fc > 10) and low spherical degree 
(I > 2) could be excited even without any enhancement 
of the iron abundance, but only in a narrow range of ef- 
fective temperature around 18 000 K. The width of the g- 
mode instability strip and the number of excited modes in- 
creases substa n tially with increasing iron abundance, but 
Ijefferv fc Said 120061) were unable to obtain unstable g- 
modes for T e g > 24 000 K. Even these modes were only ob- 
tained with iron enhanced by a factor of twenty over a base 
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Table 1. Assumed fraction of iron-grou p elements within th e 
heavy eleme nt fraction Zo: OPA L = lOrevesse Noels! Il993t) . 
OP = S92 of lSeaton et"aLNl994ly 



El 


Z 


Z e ] 


/Z 






OPAL 


OP 


Cr 


21 


0.000329 


0.000324 


Mn 


25 


0.000170 


0.00017 


Fe 


26 


0.021877 


0.02244 


Ni 


28 


0.001293 


0.00123 



metallicity Zo = 0.02 and for I = 3. This theoretical g-mode 
blue edge is some 5 000 K cooler than the hottest known 
PG1716 variables, Balloon 090100001 and HS 0702+6043 
which have T cff ~ 29 500K JOreiro et al]l200St ISchuh et all 
2006). Efforts to reduce this difference by varying the he- 
lium and iron abundances, and the overall metallicity, were 
unsuccessful. 

Subsequently the authors have considered what other 
mechanisms might contribute to the instability of g-modes 
in PG1716 stars. This paper describes the surprising results 
that the explanation lies in the atomic physics used to com- 
pute the stellar opacities (section 2) and that the theoretical 
blue-edge can be reconciled with observation by considering 
the individual roles of other iron-group elements including 
nickel (section 3). The consequences for theoretical mod- 
els of stellar interiors and the atomic processes to be found 
therein are briefly considered (section 4), before conclusions 
are drawn (section 5). 



2 STELLAR OPACITY 

The motivation for this work is the general result that the 
position of an instability strip, i.e. the T e g range within 
which a star can be expected to pulsate, is dictated by the 
temperature of the opacit y peak which provid es the driv- 
ing. In their original work. ljefferv fe Said <|2006T) considered 
variations in only helium (Y), all metals together (Zo), and 
iron (Zfc) - the latter being the dominant element within 
the iron group by a large factor. However, if radiative accel- 
eration acts on iron, it will also operate on other elements, 
particularly iron-group elements with electronic structures 
similar to that of iron. 

Our first question was, therefore, whether an increased 
contribution of elements other than iron might alter the po- 
sition (i.e. temperature) of the "Fe-bump". This was ad- 
dressed b y obtaining tables o f updated Rosseland mean 
opacities feadnell et al.l 120051) from the Opacity Project 
webserver for mixtures with X = 0.7, Zo = 0.02, but with 
elements chromium, manganese, iron and nickel individu- 
ally enhanced by a factor (/) of ten. Comparing these in 
the vicinity of the "Fe-bump" , it was evident that there are 
significant differences in the bump temperature for each ele- 
ment (Fig. . Considering the relative abundances of these 
elements (Table and that we had simply multiplied by 
ten the default fraction of each element (as prescribed by 
the respective webservers), it was also apparent that nickel 
and chromium are highly effective absorbers. Significantly, 
the bump temperature of nickel is markedly higher than that 
of iron, whilst those of manganese and chromium are sim- 



Table 2. Blue-edge (logT c g/K) of instability strip for g-modes 
in blue horizontal branch stars 
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ilar or lower. Therefore, we set out to investigate whether 
a nickel overabundance might drive g-mode pulsations at a 
higher T e g than iron. 

Since our stellar structure code is designed to work 
with OPAL opacity tables, a new set of tables was obtained 
from the OPAL webserver with nickel enhanced by a factor 
fm = 20, and other elements normal (X = 0.70, Z = 0.02). 
A pulsation stability analysis was carried out for a se- 
ries of zero-ag e horiz ontal-branch models, as described by 
IJefferv fe Said <|2006T) . These new opacities had a negligible 
effect on the stability results, being essentially the same as 
for the mixture with /m = 1 and /p e = 20. 

The reason became apparent when the OPAL and OP 
opacities were compared. The peak temperatures of the 
bumps ar e significantly different in the two sets (cf. Figs. 
1 and lO. ISeaton fe BadnelJl2004TI . The Fe-bump occurs at 
a higher temperature in OP than in OPAL. The OPAL Ni- 
bump is at practically the same temperature as the OPAL 
Fe-bump (Fig. |2J . It would be necessary to make the pulsa- 
tion calculations using OP data. 

This was achieved by obtaining single tables of Rosse- 
land mean opacity from the OP webserver for mixtures 
X = 0.70,0.35, 0.10,0.03,0.0, Z = 0.02, and reformat- 
ting th ese to match the OPAL tables using the OP utility 
opfit JSeatorill993l) . Type 2 OPAL tables were obtained 
from the OPAL webserver for the same standard mixtures, 
but including the 16 tables with enhanced C and O used 
to compute evolved stars as required by our codes. Both 
sets of tables were computed with elemental enhancements 
/p c ,Ni = (1, 1); (20, 1); (1, 20); and (10, 10). 

To use the OP tables within our codes, the OPAL tables 
were read in first. The OP tables were then read in to replace 
the corresponding tables with no enhanced C and O. This 
works because the interpolations required for modelling the 
H-rich layers of ZAHB stars make no use of the tables with 
enhanced C and O. 



3 G-MODES IN BLUE 

HORIZONTAL-BRANCH STARS 

As previously llFontaine et all2003l : I Jeffery fe Saiol2006l wc 

have assumed that the PG1716 variables are extreme hori- 
zontal branch stars and have constructed series of "zero-age" 
horizontal-branch (ZAHB) models having a helium-burning 
core with a mass of M c and with a hydrogen-rich envelope 
with mass M e ranging from 3 x 10 -6 to 0.034 M . Here 
we have considered only the case of M c — 0.486 Mq. The 
surface layers are characterized by mass fractions of hydro- 
gen X = 0.75 and base metallicity Zo = 0.02. with iron 
and nickel enhanced by the pairs of /F e ,/Ni given above. 
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Figure 1. Comparison of log(re R ) from OP for X = 0.70, Zq = 0.02 (black) with the same mixture with Cr, Mn, Fe, and Ni (left to 
right) enhanced selectively by a factor 10 (red). Curves are given for log(i?) = —2, —3, —4 and —6 (top to bottom), where R = p/Tg, p 
is mass desnity in g cm" " 3 and T 6 is 10 6 X T with T in K. 




Figure 2. Comparison of log(K R ) from OP (black) with that from OPAL (blue) for X = 0.70, Z = 0.02 (left) and with Fe (center) 
and Ni (right) enhanced by a factor 20. Otherwise as Fig. HI 



IJefferv fe Said (2006) discuss variations with X and Zq, and 
M c . In particular, varying M c has no effect on the stability 
of g-modes. They also found the only significant effect of 
evolution is that g-modes cease to be excited when the core 
becomes radiative at the end of core He-burning. Of course 
pulsation periods P increase with radius R. Since P oc R 3 ^ 2 
or surface gravity g _3//4 , normal horizontal-branch evolution 
produces a period increase of ~ 0.2 dex in all modes. 

The high-temperature end of this ZAHB sequence cor- 
responds with the locus of the short-period EC 14026 vari- 
ables. PG1716 variables are to be fou nd at lower tempera- 
tures, with 21 000 < r eff /K < 29 000 jFontaine et alJl2006t 
iRandall et al.1120 06). The lower limit probably does not rep- 
resent a formal red edge; horizontal branch stars at these 
temperatures are very scarce. 

We have tested the stability of each of our horizontal- 
branch models for non-radial p- and g-modes with spheri- 
cal degree I = 1, . . . , 4. The frequency range considered is 
0.2 < lo < 20, where ui is the angular frequency of pulsa- 
tion normalized by y/ GM/R 3 with G being the gravitational 
constant. 

Fig. |3]compares the results for OP and OPAL with nor- 
mal composition. With OP, unstable modes are obtained 
even for I = 1. In general, the number and the tempera- 
ture range of unstable modes is increased with OP opac- 
ities. Fig. [I] compares the results for OP and OPAL with 
/Fc,Ni = (10, 10). For modes with I — 3, g-modes are ex- 
cited up to T cff /K < 28 000 in OP models, but only to 
T cff /K < 25 000 in OPAL models. Although less likely to 
be observable, they are excited up to T e g /K < 29 500 for 



OP models with I — 4. In contrast, the position of the mod- 
els on the HR diagram hardly changes with the choice of 
opacity table. 

Table |2] compares the theoretical blue-edge of the g- 
mode instability strip for different enhancements of iron and 
nickel, for OP and OPAL models and for the three most 
observable modes, I = 1,2 and 3. 

It is clear that the use of OP opacities combined with 
the inclusion of excess nickel can shift the theoretical blue- 
edge of the g-mode instability strip significantly closer to 
the observed blue ed ge. This prompts us to believe that 
iFontaine et all <l2003l) correctly identified the oscillations in 
PG1716 variables as opacity-driven g-modes, and that the 
discrepancy in the predicted and observed blue edges can 
be solved by invoking more accurate atomic physics in the 
calculation of stellar opacity, as well as considering atomic 
species other than iron. 

The basic elements of this solution are that (i) the opac- 
ity peaks due to iron and nickel occur at a higher tempera- 
ture in OP than in OPAL, (ii) the OP nickel peak occurs at 
a higher temperature than the OP iron peak, (iii) at these 
temperatures nickel is a more efficient absorber than iron 
(per ion), and (iv) radiative acceleration will force it to ac- 
cumulate in the Ni-bump region in exactly the same way as 
iron. 
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Figure 3. Periods of modes due to Fe-bump instability for ZAHB stars with M c = 0.486 Mq, X = 0.75, Zq = 0.02, and M e = increasing 
from left to right. The top row shows models computed with OPAL opacities, the bottom row with OP opacities. Stable modes are 
marked as (blue) dots, unstable modes are marked as filled (red) circles. 



4 ATOMIC PHYSICS 

It should, perhaps, have been no surprise that the introduc- 
tion of OP opacities would shift the instability blue edge. A 
more modest shift in the blue edge for high-o rder g-modes 
in slow ly pulsating B stars was obtained by IPamvatnvkhl 
(1999). Since then, extensive work on the contribution of 
inner-shell electrons to the OP calculations, has brought the 
OP opac ities into "much closer agreement with thos e from 
OPAL" (|Badnell fc SeatorJl2004 iBadneU et alJl200.^ . 

Directed at obtaining good atomic data for very 
large numbers of atomic transitions, OPAL uses single- 
config uration wavefunc tions with one-electron orbitals 
^Rogers fc Ielesiasl IT992). OP originally obtained most of 
its atomic data using the R-matrix method to evalu- 
ate multi-electron wavefunctions, resulting in more accu- 
rate atomic data, but for fewer and simpler transitions 
<The Opacity Project Teamlll995lll997l) . Prompted by sub- 
stantial differences between the first OP results and OPAL, 
primarily at high temperature and densities, the introduc- 
tion of configuration-interaction codes enabled OP to in- 
clude many more bound states and radiative transition prob- 
abilities, as well as auto- ionization pro babilities and photo- 
ionization cross-sections 

For many astrophysical situations, OP and OPAL opac- 
ities app ear to be in closer agreement with one another 
(< 5%) iJSeatonlEfflflh than with observation. For exam- 
ple, helioseismology results currently require an increase of 



~ 20 % in opacity below the base of the solar convective 
zone ijBahcall et abl feoOS). While the results presented here 
do not directly suggest a solution to the solar opacity prob- 
lem, they do indicate that improvements in the treatment of 
atomic processes used in the calculation of stellar opacities 
can still have substantial consequences for calculations of 
pulsation stability. These are especially important when am- 
plified by radiatively-driven diffusion of elements into layers 
of high specific opacity, as in the outer layers of many early- 
type stars. In turn, these will impact on elemental stratifica- 
tion, on stellar radii, on boundaries of pulsational stability, 
and on pulsation periods derived for use in asteroseismology. 



5 CONCLUSION 

We have analyzed the stability of models of hot zero-age 
horizontal-branch stars to nonradial g-mode oscillations. We 
have tested and confirmed the hypothesis that elements 
other than iron can raise the temperature of the so-called 
"Fe-bump" in opacity at around 200 000 K. In doing so we 
have discovered that, within this "Fe-bump", recently up- 
dated OP opacity tables provide more opacity at signif- 
icantly higher temperatures than the OPAL tables. The 
higher temperature Fe- and Ni-bumps in the OP tables have 
been combined with a tenfold enhancement of both nickel 
and iron abundances in the driving zone. This provides suf- 
ficient opacity at the right temperatures to shift the theo- 
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retical g-mode blue-edge to within 1000 K of the observed 
blue-edge for PG1716 stars. Given the approximate nature 
of our own calculations and the likely remaining uncertain- 
ties in the atomic physics, we concl ude that, by modifyin g 
the driving mechanism identified bv lFontaine et all (2003), 
the discrepany in the predicted and observed blue edges for 
PG1716 variables can be solved. 

In reaching this conclusion, we have demonstrated that 
g-mode pulsations in PG1716 stars provide a critical test 
of stellar opacity, substantially preferring the updated OP 
calculations to OPAL. Future work should include an inves- 
tigation of the use of OP opacities in models for EC 14026 
variables and OP radiative accelerations in self-consistent 
models of chemical diffusion in extreme horizontal-branch 
stars, as well as the influence of important opacity providers 
other than iron and nickel (e.g. chromium, manganese, vana- 
dium). Detailed work for asteroseismological analyses will 
need to consider a wide range of additional physics, including 
the effects of evolution and time-dependent diffusion theory 
and core mass on predicted periods. 
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